STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS DRIVEN 

BY PURELY SPATIAL NOISE 
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Abstract. We study bilinear stochastic parabolic and elliptic PDEs driven by 
purely spatial white noise. Even the simplest equations driven by this noise often 
do not have a square-integrable solution and must be solved in special weighted 
spaces. We demonstrate that the Cameron-Martin version of the Wiener chaos 
decomposition is an effective tool to study both stationary and evolution equations 
driven by space-only noise. The paper presents results about solvability of such 
equations in weighted Wiener chaos spaces and studies the long-time behavior of 
the solutions of evolution equations with space-only noise. 
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1. Introduction 

Stochastic PDEs of the form 

ii(t, x) = Au(t, x) + Mu (t, x) ■ W (t, x) , (1.1) 

where A and M are linear partial differential operators and W (t, x) is space-time 
noise, are usually referred to as bilinear evolution SPDEsS These equations are 
of interest in various applications: nonlinear filtering for hidden Markov models [30J , 
propagation of magnetic field in random flow [2] stochastic transport [SI [20]), porous 
media [3] and others. The theory and the applications of bilinear SPDEs have been 
actively investigated for a few decades now; see, for example, [U [13], [231 123 (211 [29] 
etc.]. 

In contrast, very little is known about bilinear parabolic and elliptic equations driven 
by purely spatial Gaussian white noise W (x) . Important examples of these equations 
include: 
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Bilinear SPDEs differ from linear by the term including multiplicative noise. Bilinear SPDEs 
are technically more difficult then linear. On the other hand, multiplicative models preserve many 
features of the unperturbed equation, such as positivity of the solution and conservation of mass, 
and are often more "physical" . 
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1. Heat equation with random potential modeled by spatial white noise: 

u(t, x) = Au(t, x)+u (t, x)oW (x) , (1.2) 

where o denotes the Wick product, which, in this case, coincides with the Skorokhod 
integral in the sense of Malliavin calculus. A surprising discovery made in [9J was 
that the spatial regularity of the solution of equation (jl.2p is better than in the case 
of similar equation driven by the space-time white noise. 

2. Stochastic Poisson equations in random medium: 

V(A £ (x)oVu(x)) = f(x), (1.3) 

where A £ (x) := (a(x) + eW (x)) , a(x) is a deterministic positive-definite 
matrix, and e is a positive number. (Note that a(x) o Vii (x) = a(x)Vu (x)). 

3. Heat equation in random medium: 

v (t, x) = V {A £ (x) o Vv (t, x)) + g(t, x). (1.4) 

Note that the matrix A e in equations (jl.3p and (jl.4p is not necessarily positive definite; 
only its expectation a (x) is. 

Equations (11.31) and (11.41) are random perturbation of the deterministic Poisson and 
Heat equations. An important feature of these type of perturbation is that the result- 
ing equations are unbiased in that they preserve the mean dynamics: the functions 
uo (x) := Eit (x) and vq := Et> (t, x) solve the deterministic Poisson equation, 

V(a(x)Vw (t,x)) = Ef (x) 
and the deterministic Heat equation 

Vq (t, x) — V (a (x) Vq (t, x)) + ~Kg(t, x), 

respectively. 

The objective of this paper is to develop a systematic approach to bilinear SPDEs 
driven by purely spatial Gaussian noise. More specifically, we will investigate bilinear 
parabolic equations 

= Av (t, x) + Mv(t, x) o W(x) - f(x) (1.5) 

and elliptic equations 

Au(i) + Mu(i)^(i) = /(i), (1.6) 
for a wide range of operators A and M. 

Purely spatial white noise is an important type of stationary perturbations. However, 
except for elliptic equations with additive random forcing [5l [221 [26] . SPDEs driven 
by spatial noise have not been investigated nearly as extensively as those driven by 
strictly temporal or space-time noise. 

In the case of spatial white noise, there is no natural and convenient filtration, es- 
pecially in the dimension d > 2. Therefore, it makes sense to consider anticipative 
solutions. This rules out Ito calculus and makes it necessary to rely on more nuanced 
techniques provided by Malliavin calculus. 
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In this paper, we deal with broad classes of operators A and M that were investigated 
previously for causal (nonanticipating) solutions of the equation (II. ip driven by space- 
time white noise. 

The notion of ellipticity for SPDEs is more restrictive then in deterministic theory. 
Traditionally, nonanticipating solutions of equation (11.11) were studied under the fol- 
lowing assumptions: 

(i) The operator A — |MM* is "elliptic" (possibly degenerate coercive) operator. 

Of course, this assumption does not hold for equation (jl.4p and other equations in 
which the operators A and M have the same order. Therefore, it is important to 
study equations (11.61) and (II. 5p under weaker assumptions, for example 

(ii) The operator A is coercive and ord (M) < ord (A) . 

In 1981, it was shown by Krylov and Rozovskii [T3] that, unless assumption (i) holds, 
equation (11.11) has no solutions in the space L2{VL\X) of square integrable (in prob- 
ability) solutions in any reasonable functional space X. The same effect holds for 
bilinear SPDEs driven by space only white noise. 

Numerous attempts to investigate solutions of stochastic PDEs violating the sto- 
chastic ellipticity conditions were made since then. In particular it was shown in 
[121 [25] that if the operator A is coercive ("elliptic") and ord (M) is strictly less 
then ord (A) then there there exists a unique generalized (Wiener Chaos) nonantici- 
pative solution of equation (II. ip . This generalized solution is a formal Wiener Chaos 
series u = J2\ a \<oo u ^c^ where {£ a }| a | <00 is the Cameron-Martin orthonormal basis 
in the space /^(fi). Regularity of this solution is determined by system of positive 
weights {r a }\ a \ <00 and a function space X such that 



The stochastic Fourier coefficients u a satisfy a lower-triangular system of deterministic 
PDEs. This system, called propagator, is uniquely determined by the underlying 
equation (11. ip . 

Stochastic spaces equipped with the norms similar to (11. 7B have been known for quite 
some time; see, e.g. [ITJ [T2J ES]- For historical remarks regarding other types of 
generalized solutions and applications to SPDEs see the review paper [18] and the 
references therein. 

The Wiener Chaos is a bona fide generalization of the classical Ito solution: if exists, 
a non-antisipating square integrable Ito solution coincides with the Wiener Chaos 
solution. 

In this paper, we establish existence and uniqueness of Wiener Chaos solutions for 
stationary (elliptic) equations of the type ( 11.61) and evolution (parabolic) equations 
of the type ( 11. 5ft . These results are proved under assumption (ii) that allows us in 
particular to deal with equations like ( 11. 3D and ( 11. 4ft . In many cases we are able to 
find optimal or near-optimal systems of weights {?"a}| a |<oo that guarantee (I1.7P 




(1.7) 



|«|<oo 
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Finally, we establish the convergence, as t — > +00, of the solution of the evolution 
equation to the solution of the related stationary equation. 

The structure of the paper is as follows. Section [2] reviews the definition of the Sko- 
rokhod integral in the framework of the Malliavin calculus and shows how to compute 
the integral using Wiener chaos. Sections [3] and H] deal with existence and uniqueness 
of solutions to abstract evolution and stationary equations, respectively, driven by a 
general (not necessarily white) spatial Gaussian noise; Section @] describes also the 
limiting behavior of the solution of the evolution equation; Section illustrates the 
general results for bilinear SPDEs driven by purely spatial white noise. 

2. Weighted Wiener Chaos and Malliavin Calculus 

Let F = (fi, J 7 , P) be a complete probability space, and U, a real separable Hilbert 
space with inner product (-, -)u- On F, consider a zero-mean Gaussian family 

W = {W(h) } heh(} 

so that 

E(w(/n) W(h 2 ^ = (h u h 2 )u. 

It suffice, for our purposes, to assume that T is the cx-algebra generated by W. Given 
a real separable Hilbert space X, we denote by L 2 (W; X) the Hilbert space of square- 
integrable T- measurable X- valued random elements /. In particular, 

(f,9)Uw;X)-=m9)x- 

When X = R, we write L 2 (F) instead of L 2 (¥; R). 
Definition 2.1. A formal series 

W = Y^W(u k )u k , (2.1) 

k 

where {u k , k > 1} is a complete orthonormal basis in U, is called (Gaussian) white 
noise on 14. 

The white noise on U = L 2 (G) , where G is a domain in R d , is usually referred to 
as a spatial or space white noise (on L 2 (G)). The space white noise is of central 
importance for this paper. 

Below, we will introduce a class of spaces that are convenient for treating nonlinear 
functionals of white noise, in particular, solutions of SPDEs driven by white noise. 

Given an orthonormal basis it = {u k ,k > 1} in U, define a collection k > 1} 
of independent standard Gaussian random variables so that = W{\ik). Denote 
by J the collection of multi-indices a with a = (a±,a 2 , . . .) so that each a k is a 
non-negative integer and \a\ := J2k>i a k < 00 • F° r oz, (3 G J, we define 

a + (3 = (cti + Pi,a 2 + (h, ■ ■ .), a! = J^[a fc !. 

fe>i 
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By (0) we denote the multi-index with all zeroes. By E{ we denote the multi-index a 
with aj = 1 and otj = for j ^ i. With this notation, nei is the multi-index a with 
a, = n and aj = for j ^ i. The following inequality holds (see Appendix for the 
proof): 

\a\\ < a!(2N) 2Q , where (2N) 2 ° = JJ(2A;) 2cifc . (2.2) 

k>l 

Define the collection of random variables H = a G JT} as follows: 

«- = Il(^#). < 2 - 3 > 



where 

H n (x) = (-l) n e* 2/2 ^e-* 2 / 2 (2.4) 
is Hermite polynomial of order n. 

Theorem 2.2 (Cameron and Martin [1]). The collection 5 = {£ a , a G J7"} an 
orthonormal basis in L 2 (F): if i] E L2OF) and = E(r/^ a ), i/ien 77 = J] Q ^ a and 

Expansions with respect to the Cameron-Martin basis H is usually referred to as 
Wiener Chaos. Next, we introduce a modification of the Wiener Chaos expansion 
which will be called weighted Wiener Chaos. 

Let 7Z be a bounded linear operator on L 2 (F) defined by lZ£, a = r a £ a for every a <E J , 
where the weights {r a , ol G J} are positive numbers. By Theorem 12.21 71 is bounded 
if and only if the weights r a are uniformly bounded from above: r a < C for all a G J~, 
with C independent of a. The inverse operator Tlr 1 is defined by lZ~ l S, a = r" 1 ^. 

We extend 1Z to an operator on L 2 (F;X) by defining IZf as the unique element of 
L 2 (F;X) so that, for all g G £ 2 (F;X), 

E(ft/,0) x = £r a E((/,0) x £«). 
Denote by 7?X 2 (F; X) the closure of L 2 (F; X) with respect to the norm 

ll/|lwi 2 (F;Jf) := II^/IIl2(F;X)- 

Then the elements of 7?X 2 (F;X) can be identified with a formal series Yla^j fa£,a, 
where f a eX and Y, a ej \\f»\\x r l < °°. 

We define the space 7^ _1 L 2 (F;X) as the dual of 7£L 2 (F;X) relative to the inner 
product in the space L 2 (M; X) : 

ft-%(F;X) = {ge L 2 (¥;X) : TT x g G L 2 (¥;X)} . 

For / G 7?X 2 (F; X) and g G 7?.^ 1 L 2 (F) we define the scalar product 



«/, g)):= ¥,({Kf){n- l g)) el 



(2.5) 
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In what follows, we will identify the operator 1Z with the corresponding collection 
(r a , a G J). Note that if u G 1ZiL 2 (¥; X) and v G TZ 2 L 2 (¥;X), then both u and v 
belong to TZL 2 (¥; X), where r a = min(r ljQ ,, r 2jQ ,). As usual, the argument X will be 
omitted if X = JR. 

Important particular cases of the space 1ZL 2 (¥;X) correspond to the following 
weights: 

(1) 

oc 

*:=i 

where {qk, k > 1} is a non-increasing sequence of positive numbers with q x < 1 
(see USES]); 

(2) 

r 2 a = (a\) p (2NY a , p < 0, £ < 0, where (2N) to = JJ(2A;) tofc . (2.6) 

k>l 

This set of weights defines Kondratiev's spaces (S) P ^(X). 

Now we will sketch the basics of Malliavin calculus on 7^X2 (F; X). 

Denote by D the Malliavin derivative on L 2 (¥) (see e.g. [24j ) . In particular, if 
F : W N — > R is a smooth function and hi G W, i = 1, . . . N, then 

DF(W(/n), . . . t^(M) = ■ ■ ■ , W(h N ))hi G L 2 (F; W). (2.7) 

i=l 1 

It is known [21] that the domain © 1,2 (F) of the operator D is a dense linear subspace 
of /-2(F). 

The adjoint of the Malliavin derivative on L 2 (W) is the Ito-Skorokhod integral and is 
traditionally denoted by 5 |24j. We will keep this notation for the extension of this 
operator to TZL 2 (¥; X <g> U). 

For / G TZL 2 (¥; X® U), we define 5(f) as the unique element of 7ZL 2 (¥; X) with the 
property 

((5(f), V }}=E(1lf,'Jl- l Vip)u (2.8) 
for every ip satisfying ip G 7£. _1 L 2 (F) and Dip G TZ~ 1 L 2 (¥;U). 

Next, we derive the expressions for the Malliavin derivative D and its adjoint 5 in 
the basis 5. To begin, we compute D(£ a ). 

Proposition 2.3. For each a G 3 , we have 

D(£ a ) = £Va^ a _ efc u*. (2-9) 
fc>i 



Proof. The result follows by direct computation using the property (I2.7P of the 
Malliavin derivative and the relation H' n (x) = nH n _i(x) for the Hermite polynomials 

(cf. ES|). ' □ 
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Obviously, the set J is not invariant with respect to substraction. In particular, the 
expression a — e k is undefined if a k = 0. In (12.91) and everywhere below in this paper 
where undefined expressions of this type appear, we use the following convention: if 
a k = 0, then y/a^^ a -e k = 0. 

Proposition 2.4. For £ a G H, h G X, and u k G it, we have 



5(£ah®U k ) = hV^k + l£a+e k . (2.10) 

Proof. It is enough to verify (12.81) with / = h <E> u k £ a and ip = £p, where h G X. By 

(ED, 



VOfcTl fr, if a = /? - ejfc, 
0, ifa^P- e k . 



In other words, 

E(£ Q /i <g> u fe , D^)w = ftE(Va, + lU&) 
for all f3 <E J. □ 

Remark 2.5. T/ie operator <5D zs linear and unbounded on L 2 (F); i£ follows from 
Propositions \2. 31 and IL^ i/iai £/ie random variables £ a are eigenf unctions of this op- 
erator: 



6(D(£ a )) = \a\U (2.11) 

To give an alternative characterization of the operator 5, we define a new operation 
on the elements of S. 

Definition 2.6. For £ a; ^ /rom H, define the Wick product 



^o^:=W(^p)w (2-12) 

In particular, taking in (12.61) a = fee, and (3 = ne^ and using (I2.3p . we get 

# fc (&)o# n (&) =H k+n {ii). (2.13) 

By linearity, we define the Wick product for] for / G 7£L 2 (F; X) and 77 G 7^L 2 (F): 
if / = Ysaaj fc&m fa e -X", and 77 = E q6 j ?7c*£a, e R, then 

Proposition 2.7. If f G 7?X 2 (F;X) and 77 G 7?X 2 (F), £/ien f o r] is an element of 
TZL 2 (¥; X) /or a suitable operator 71. 

Proof. It follows from (12. 6p that / o 77 = Xlaej ^a£<* an d 
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Therefore, each F a X is an element of X, because, for every a G J , there are only 
finitely many multi-indices /3, 7 satisfying (3 + 7 = a. It is known [TBI Proposition 
7.1] that 

^(2N)" Q < 00 if and only if q < -1. (2.14) 

Therefore, f orj G 1ZL 2 (¥; X), where the operator TZ can be defined using the weights 
rl = (2N)- 2 7(1 + ||F Q ||^). ' □ 

An immediate consequence of Proposition 12.41 and Definition 12.61 is the following 
identity: 

S{£ a h <g> u fc ) = h£ a o £ fc , he X. (2.15) 
Below we summarize the properties of the operator S. 

Theorem 2.8. If f is an element of TZL 2 (¥; X <g> hi) so that f = J2 k>1 fk <£> ^k, with 
fk = J2aej fk,a£a e TIL 2 (¥; X) , then 

*(/) = ( 2 ' 16 ) 

k>l 

and 

W)) a = J2V^fk,a- ek . (2.17) 

fc>l 



Proof. By linearity and (I2.15p . 

s (f) = ^2^2S{£ a fk,a®u k ) = ^2^2f k ,a£ a o£k = ^2fk<>£k, 

k>\ a&J k>l aej k>l 

which is (12161) . On the other hand, by (l2~T0i) . 

8(f) = ^2^2 fk,aV(XkTl£a+e k = fk,a-e k V^k Ca, 

k>l a£j k>l aej 

and OTTj) follows. □ 

Remark 2.9. zs noi difficult to show that the operator S can be considered as an 
extension of the Skorokhod integral to the weighted spaces TZL 2 (W; X ® 14). 



One way to describe a multi-index a with \a\ = n > is by its characteristic set K a , 
that is, an ordered n-tuple K a = {hi, . . . , k n }, where k\ < k 2 < . . . < k n characterize 
the locations and the values of the non-zero elements of a. More precisely, k\ is 
the index of the first non-zero element of a, followed by max (0, a kl — 1) of entries 
with the same value. The next entry after that is the index of the second non-zero 
element of a, followed by max (0, a k2 — 1) of entries with the same value, and so on. 
For example, if n = 7 and a = (1, 0, 2, 0, 0, 1, 0, 3, 0, . . .), then the non-zero elements 
of a are ot\ = 1, 03 = 2, a§ — 1, as — 3. As a result, K a = {1,3, 3, 6, 8, 8, 8}, that is, 
k\ = 1, k 2 = k 3 = 3, &4 = 6, k§ = k% = k 7 = 8. 
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Using the notion of the characteristic set, we now state the following analog of the 
well-known result of Ito [TU] connecting multiple Wiener integrals and Hermite poly- 
nomials. 

Proposition 2.10. Let a G J be a multi-index with \a\ — n > 1 and characteristic 
set K a = {k%, . . . , k n }. Then 

£ Q = 61 062 (2.18) 
Va! 

Proof. This follows from (12.31) and (I2.13p . because by f)2.13p . for every % and k, 

H k {£,i) = & o • • • o & . 

A; times 

□ 



3. Evolution Equations Driven by White Noise 

3.1. The setting. In this section we study anticipating solutions of stochastic evo- 
lution equations driven by Gaussian white noise on a Hilbert space 14. 

Definition 3.1. The triple (V, H, V) of Hilbert spaces is called normal if and only 
if 

(1) V ■=— > H <^-> V and both embeddings V ^ H and H •— > V are dense and 
continuous; 

(2) The space V is the dual of V relative to the inner product in H; 

(3) There exists a constant C > so that \(h,v)n\ < C||u||y||/i||y/ for all v G V 
and h G H . 

For example, the Sobolev spaces (H e 2 + "'(R d ), #|(M d ), ^ 7 (M d )), 7 > 0, £ G R, form 
a normal triple. 

Denote by (v',v), v' G V, f G V, the duality between V and V relative to the inner 
product in H. The properties of the normal triple imply that \ (v',v)\ < C\\v\\v\\v'\\v>, 
and, if v ' G H and v G V, then = (i/,t>)#. 

We will also use the following notation: 

V = L 2 ((0,T);\/), n = L 2 ((0,T);H), V = L 2 ((0,T); V). (3.1) 

Given a normal triple (V, if, V), let A : V — »• V and M : V -> V" <g> W be bounded 
linear operators. 

Definition 3.2. TTie solution of the stochastic evolution equation 

u = Au + f + 5(Mu), < t < T, (3.2) 

wzt/i / G 7?X 2 (F; V) and u\ t =o = u G 7?X 2 (F; if), is a process u G 7?X 2 (F; V) so £/ia£, 
for every ip satisfying tp G 7?. _1 L 2 (F) and D(/? G 7?. _1 L 2 (F; W) ; £/ie equality 

((u(t),<p)) = ((uo,<p)) + [\au(s) + f(s) + 8(Mu)(s),cp))ds (3.3) 

Jo 
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holds in V' ; see ( QL 51) for the definition of ((•, •, )). 

Remark 3.3. (a) The solutions described by Definitions \ 3.2\ and \4-l\ belong to the 
class of "variational solutions", which is quite typical for partial differential equations 

(see nam EM, etc -) 

(b) Since ((u(t),ip)) G V and ((u(t),ip)) t G V, by the standard embedding theorem (see 
e.g. [TO Section 1-2.2}) there exists a version of ((u(t),<p)) G C ([0,T]; H). Clearly, 
one could also select a version of u{t) such that ((u(t),(p)) G C ([0, T]; H) . In the 
future, we will consider only this version of the solution. By doing this we ensure that 
formula \3. 3\) which is understood as an equality in V yields u\ t =o = u$ G 7^X 2 (F; H). 

Remark 3.4. To simplify the notations and the overall presentation, we assume that 
A and M do not depend on time, even though many of the results in this paper can 
easily be extended to time- dependent operators. 

Fix an orthonormal basis il in U. Then, for every v G V, there exists a collection 

Vk G V, k > 1, so that 

M.v = 2J Vf. <S> u k . 

We therefore define the operators : V — > V by setting M^f = and write 

Mv = y2(M k v) ®u fe . 

k>l 

By (12.161) . equation f)3.2p becomes 

u(t) = Au(t) + f(t)+Mu(t)oW, (3.4) 

where 

MvoW:=^2{M k v)o^ k . (3.5) 

k>\ 



3.2. Equivalence Theorem. In this section we investigate stochastic Fourier rep- 
resentation of equation (13 .4p . 

Recall that every process u = u(t) from 7^X2 (F; V) is represented by a formal series 
u(t) = J2 aGj u a (t)£ a , with u a G V and 

^2 r l IKHv < 00. (3.6) 

a 

Theorem 3.5. Let u = Y2 a ej Ua £ a ^ e an e/emeni ofTZL>2(F; V). The process u is a 
solution of equation / tff.jgj) if and only if the functions u a have the following properties: 

(1) every u a is an element of C ([0, T}; H)) 

(2) the system of equalities 

u a (t) = u , a + A.u a (s) + f a (s) + ^ x /q^M fc tt Qi - £fc (g) ) ds (3.7) 
70 V k>i J 

holds in V for all t G [0, T] and a G J . 
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Proof. Let u be a solution of (13.21) in 7^L 2 (F; V). Taking <p = £ a in (13.31) and using 
relation fl2TT7j) . we obtain equation (JXTD- By Remark 6 Vf|C ([0,T];H)) . 

Conversely, let {u a , a G J7"} be a collection of functions from V f] C ([0, T]; if)) satis- 
fying (I3.6P and (13.71) . Set u (t) := J2aej u ait)^ a . Then, by Theorem I2.8[ equation 
(13.71) yields that, for every a G J , 

= KU+ A(A«( S )+/( S )+«J(Mu)(s),e a ))ds. 







By continuity, we conclude that for any <p ElZ 1 L 2 (¥) such that D<^ G 7?. 1 L 2 (F;W), 
equality 

««(*),¥>» = + [ ((Au(s)+f{s) + 6(Mu)(s),<p))da 

Jo 

holds in V. By Remark E£2 y?)) G C ([0, T]; if) . □ 

This simple but very helpful result establishes the equivalence of the "physical" (13.41) 
and the (stochastic) Fourier (13.71) forms of equation (13.21) . System of equations (13.71) 
is often referred in the literature as the propagator of equation (13.41) . Note that the 
propagator is lower-triangular and can be solved by induction on \a\. 

3.3. Existence and uniqueness. Below, we will present several results on existence 
and uniqueness of evolution equations driven by Gaussian white noise. 

Before proceeding with general existence-uniqueness problems, we will introduce two 
simple examples that indicate the limits of the "quality" of solutions of bi-linear 
SPDEs driven by general Gaussian white noise. 

Example 3.6. Consider equation 

t 



u(t) = <f>+ / (bu(s)o£- Xu(s))ds, (3.1 







where <f), A are real numbers, b is a complex number, and £ is a standard Gauss- 
ian random variable. In other words £ is Gaussian white noise on U = R. With 
only one Gaussian random variable the set J becomes {0, 1,2,.. .} so that u(t) = 
^2 n>0 U( n )(t)H n (^)/ \fn\, where H n is Hermite polynomial of order n According 

to fry, 



U(n)(t) = 4>I(n=0) ~ \ \U( n )(s)ds+ / byfnU(n-l)(s)ds. 

Jo Jo 

It follows that M(o)(^) = 4>e~ xt and then, by induction, u^ n )(t) = <p — i=e~ xt ■ As a 
result, 



u(t)=e-*(<f ) + J2^f H ^ 



^ e -Xt+(bt^\b\ 2 t 2 /2)^ 



Obviously, the solution of the equation is square integrable on any fixed time interval. 
However, as the next example indicates, the solutions of SPDEs driven by stationary 
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noise are much more intricate then the non-anticipating, or adapted, solutions of 
SPDEs driven by space-time white noise. 

Example 3.7. With £ as in the previous examples, consider a partial differential 
equation 

u t (t, x) = au xx (t, x) + (/3u (t, x) + <ju x (t, x)) o £, t > 0, iGR, (3.9) 

with some initial condition Uq G L 2 (M). By taking the Fourier transform and using 
the results of Example \3.6\ with <fi = Uo(y), A = —ay 2 , b = (3 + \f— lye, we find 

u t (t, y) = -y 2 au (t) + ((3 + \f-\yo) u (t, y) o f ; 

u(t, y) = u {y) exp {-tay 2 + (a 2 y 2 - /3 2 ) t 2 /2 + V^T(3ayt 2 + {V^lay + 0) t£) . 

If a = 0, i.e. the "diffusion" operator in equation $3. 9\) is of order zero, then the 
solution belongs to L 2 (F; ^(K.)) for all t. However, if a > 0, then the solution u(t, •) 
will, in general, belong to L 2 (F; L 2 (M)) only for t < 2a/a 2 . This blow-up in finite 
time is in sharp contrast with the solution of the equation 

■u t = au xx + au x o w, (3.10) 

driven by the standard one- dimensional white noise w (t) = d t W (t), where W (t) is 
the one- dimensional Brownian motion; a more familiar way of writing A3.10\) is in 
the ltd form 

du = au xx dt + au x dW(t). (3- 11) 

It is well known (see, for example, [28\) that the solution of 113.11]) belongs to 
L 2 (F; L 2 (IR)) for every t > as long as Uq G L 2 (M.) and 

a-CT 2 /2>0. (3.12) 

The existence of a square integrable (global) solution of an Ito's SPDE with square 
integrable initial condition hinges on the parabolic condition which in the case of 
equation f)3.10p is given by f)3.12p . Example 13.71 shows that this condition is not in 
any way sufficient for SPDEs involving a Skorokhod-type integral. The next theorem 
provides sufficient conditions for the existence and uniqueness of a solution to equation 
(13.41) in the space 7£L 2 (F; V), which appears to be a reasonable extension of the class 
of square integrable solutions. 

Firstly, we introduce an additional assumption on the operator A that will be used 
throughout this section: 

(A): For every Uq G H and F G V' := L 2 ((0,T); V), there exists a function U G V 
that solve the deterministic equation 

d t U(t) = AU(t) + F(t), U (0) = U Q , (3.13) 

and there exists a constant C = C (A,T) so that 

||[/|| v <CXA,T)(||[/ ||tf + ||F|| v /). (3.14) 



STOCHASTIC EQUATIONS WITH SPATIAL NOISE 13 

Remark 3.8. Assumption (A) implies that a solution of equation $3.13\) is unique 
and belongs to C ((0, T);if) (cf. Remark \3.$) . The assumption also implies that the 
operator A generates a semi-group $ = t > 0, and, for every v G V, 

2 

dt <Cl\\v\\ 2 v , (3.15) 

v 



$ t _ s M fc f (s) ds 
with numbers Cf. independent of v. 



Remark 3.9. There are various types of assumptions on the operator A that yield 
the statement of the assumption (A). In particular, (A) holds if the operator A is 
coercive in (V, H, V): 

(Av,v) + ^vWl < C\\v\\ 2 H 
for every v G V , where 7 > and C G M are both independent of v. 

Theorem 3.10. Assume(A). Consider equation ( [ff.^l ) in which uq G 71L 2 (W; H), 
f G TZL 2 (¥] V) for some operator TZ, and each is a bounded linear operator from 
V toV. 

Then there exist an operator 1Z and a unique solution u G TZL 2 {¥\ V) of jS.J^ . 

Proof. By Theorem l3.5[ it suffices to prove that the propagator (I3.7I) has a unique so- 
lution {u a (t)) ae j such that for each a, u a G V f] C ([0, T] ; H) and u := J2aej u a£,a £ 
1ZL 2 (¥;V). 

For a = (0), that is, when \a\ = 0, equation (I3.7I) reduces to 

u (0) (t) = w ,(o) + / (A«(o) + /(o)) (s)ds. 
Jo 

By (A), this equation has a unique solution and 

||w(o)||v < C(A,T) (||w ,(o)||h + ||/(o)l|v) • 
Using assumption (A), it follows by induction on \a\ that, for every a G J , equation 

d t u a (t) = Au a (t) + f a (t) + y^ jV ^M k u a - ek (t) , u a (0) = u 0: a (3.16) 

k>l 

has a unique solution in V P) C ([0, T] ; H) . Moreover, by (13. 14ft . 
IKHv < C(A,M,T) || 

\f(^k\\ u a-e k ||V 

V k>i 

Since only finitely many of a k are different from 0, we conclude that ||u Q ||v < 00 for 
all a G J. 

Define the operator 1Z on L 2 (W) using the weights 

(2N)" KQ 



1 + \\u a \ 



where k > 1/2 (cf. (12.61) ). Then u{t) := J2aej u aif)ia is a solution of (13.41) and, by 
fl2~14j) . belongs to 1ZL 2 (¥; V). ^ ^ □ 
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While Theorem 13.101 establishes that under very broad assumptions one can find an 
operator 1Z such that equation (13.41) has a unique solution in 1ZL 2 (F;V), the choice 
of the operator 1Z is not sufficiently explicit (because of the presence of ||w Q ||v) an d 
is not necessarily optimal. 

Consider equation H3.4[) with non-random / and u . In this situation, it is possible 
to find more constructive expression for r a and to derive explicit formulas, both for 
IZu and for each individual u a . 

Theorem 3.11. If Uq and f are non-random, then the following holds: 

(1) the coefficient u a , corresponding to the multi-index a with \a\ = n > 1 and 
characteristic set K a = {k\, . . . , k n }, is given by 

_^ rt rs„ r-S2 

U *( t ) = -F=i^2 / •••/ ®t- Sn Mk a(n) ---®s 2 -s 1 M kiT(1) u ( o)(s 1 )ds l ...ds n , 
v a! agPn Jo Jo Jo 

(3.17) 



where 

• V n is the permutation group of the set (1, 

• $t is the semigroup generated by A; 

• u (0 )(t) = <& t uo + f*^t-sf(s)ds. 

(2) the weights r a can be taken in the form 

a 00 



where q a = JJC (3.18) 



k=l 



and the numbers q^, k > 1, are chosen so that Ylk>i1k^k < 1; with C k from 
(3) With q k and r a from $3.180 . we have 



J2 1 a Mt)^ = ■■■ $ t _ Sn 5(M$ Sn _ Sn _ 1 5(. . . 8(Mu {0) )) . . .)ds x . . . ds n -xds n , 

\ a \=n J ° J ° J ° 

_ (3.19) 
where M = (^Mi, q^M^, ■ ■ ■), o,nd 

TZu{t) = U( )(*) 

°° 1 /"* f Sn f S2 

+ J2 o / / • • • / ^sJ(^s n -s^A- • • S(Mu (si))) . . .)dsi . . . ds n ^ds n . 
~ 2 n Vnl Jo Jo Jo 

(3.20) 



Proof. If uo and / are deterministic, then equation (13.71) becomes 

U(p)(t)=u + / Au(o)(s)ds+ / f(s)ds, \a\ = 0; (3.21) 
Jo Jo 

u a (t) = / Au a {s)ds + 2_]^/ c ^k / M.kU a - £k (s)ds, \a\ > 0. (3.22) 
Jo k>l Jo 
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Define Then W( ) = «(o) and, for \a\ > 0, (13.221) implies 



u, 



t (t) = / Au a (s)ds + / a k M k u a _ £k (s)ds 
Jo k>1 Jo 



or 

u, 



{t) = ^a k [ ^ s M k u a _ £k {s)ds = [ $t- s M k u a _ Ek (s)ds. 

k>l J ° keK a J ° 



By induction on n 



'<*(*)= E/ / •••/ ®t- Sn M K{n y--$ S2 „ sl M ki7m u (0) ds 1 ...ds n , 

JO JO JO 



and (I3.17P follows. 

Since (13.201) follows directly from (13.191) . it remains to establish (I3.19p . To this end, 
define 

Unit) = <l a u*(t)Z<x, n>0. 

\a\=n 

Let us first show that, for each n > 1, U n G L 2 (W; V). Indeed, for a = (0), u a (0) = u , 
f a = f and 

u (0 )(t) = $t« + / $ t -,f(s)ds. 
Jo 

By dSHD, we have 

|| U(0 )||v<C7(A,r)(||Mo||H+||/||vO- (3-23) 
When | a \ > 1, / a = and the solution of (I3.22p is given by 

u a {t) = J2V®k~ [ ®t-sM k u a - Ek (s)ds. (3.24) 
fc >i Jo 

By (T3TT7D . together with (13TT4D . f[3T23|) . and fl3TT5|) . we have 

\\uX<C\A,T)^(\\u o r H +\\f\\l,) (3-25) 



By the multinomial formula, 



a! 

fc>i 



<fc>l / |a|=n \ ' fe>l 



Then 



E? 2a iKHv<c 2 (A,T)(|Klll, + ||/|iy^ E fen^H 

|a|=n |a|=n V ' k>l / 

= c 2 (A,T)(\\u \\ 2 H + ii/nym (e^ 2 ) <00 ' 

\fe>i / 
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because of the selection of q k , and so U n G L 2 (¥; V). Moreover, if the weights r a are 
defined by (ETTB1) . then 

E^IKHv = J2 J2 rl\\u a \\l < C 2 (A,T)(|M^ + \\fWl,) J2 (j^Clql) < oc 

a&J n>0 \a\=n n>l \k>l J 

because of the assumption J2 k>1 C k q k < 1. 

Next, the definition of U n (t) and (13.241) imply that (13 . 191) is equivalent to 

U n {t) = [ ^ s S(MU n ^(s))ds, n>l. (3.27) 
Jo 

Accordingly, we will prove (I3.27p . For n = 1, we have 

U i( s ) =J2l kU ck( t )€k = J2 [ q k ®t-sM k u i0) € k dt = [ $t- s 5(MU {s))ds, 
fc>i fc>i Jo Jo 

where the last equality follows from (I2.16p . More generally, for n > 1 we have by 
definition of U n that 

_ J q a u a (t), if |a| = n, 
I U, otherwise. 

From the equation 



q a u a (t) = / Aq a u a (s)ds + / Qky^M k q a £k u a - £k (s)ds 
Jo k>1 Jo 



we find 

(U n (t)) a 



{^2Vak~qk ® t - s M k q a £k u a - Ek (s)ds, if|a|=n, 
0, otherwise. 

= ^2V®k [ $<- s M fc ([4_i(s)) Q _ £fc ds, 
k>i Jo 

and then ( 13.271) follows from (12. 1TI) . Theorem 14.51 is proved. □ 

Formula (13.191) is similar to the multiple Wiener integral representation of the solution 
of a stochastic parabolic equation driven by the Wiener process; see [T9l Theorem 3.8]. 



Example 3.12. Consider the equation 

u(t,x) = u (x) + / u xx (s,x)ds + y~] I cr k u xx (s, x) o£ k ds. (3.28) 
Jo — 1 Jo 



k>l 



With no loss of generality assume that o k ^ for all k. Standard properties of the 
heat kernel imply assumption (A) and inequality $3.15\) with C k = a\. Then the 
conclusions of Theorem \3 . 1 1\ hold, and we can take q\ = k~ 2 A~ k {I + a 2 k )- k . Note that 
Theorem \3 . 1 1\ covers equation A3. 28\) with no restrictions on the numbers o~ k . 
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In the existing literature on the subject, equations of the type (13.41) are considered 
only under the assumption 

(H): each is a bounded linear operators from V to H. 

Obviously this assumption rules out equation f!3.28j) but still covers equation (13.91) . 

Of course, Theorem 13.111 does not rule out a possibility of a better-behaving solution 
under additional assumptions on the operators Indeed, it was shown in [T8] that 
if (H) is assumed and the space-only Gaussian noise in equation ( 13. 4ft is replaced by 
the space-time white noise, then a more delicate analysis of equation (13.4D is possible. 
In particular, the solution can belong to a much smaller Wiener chaos space even if 
u and / are not deterministic. 

If the operators are bounded in H (see e.g. equation (13. 9p with o = 0), then, as 
the following theorem shows, the solutions can be square integrable (cf. [9]). 

Theorem 3.13. Assume that the operator A satisfies 

(Av,v) + k\\v \\v < C A \\v\\% (3.29) 

for every v G V , with k > 0, C A G R independent of v, and assume that each is 
a bounded operator on H so that ||Mfc| < Ck and 

C M :=J2 c l < 00 - ( 3 - 30 ) 

k>l 

If f G V and uq G H are non-random, then there exists a unique solution u of flff.^l ) 
so that u(t) G L 2 (¥; H) for every t and 

n<t)f H <C{C A ,C M ,K,t) (^V(s)ilv^+NII#) • (3.31) 

Proof. Existence and uniqueness of the solution follow from Theorem 13.101 and 
Remark 13.91 an d it remains to establish (I3.3ip . 



i H r* 

-= V / $ t _ s MkU a - £k (s)ds, (3.32) 
la\ Jo 



It follows from (13.71) that 

i 

>/',,! 

k&K a 

where $ is the semi-group generated by A and K a is the characteristic set of a. 
Assumption (13.291) implies that < e pt for some p G R. A straightforward 

calculation using relation f!3.32j) and induction on |a| shows that 

'~<>(0)\\H, 



\u a (t)\\ H <<^—=\\u, 0) \\ H , (3.33) 



'Oil 

-t 



where c a = Ylk c k k an< ^ M (o)(^) = ^t M o + f ^t-sf(s)ds. Assumption (13.29)) implies 
that ||u(o)||jj < C(C A ,K,t)(tt ||/(s)||v^ + ||«o|||)- To establish QSM), it remans 



to observe that 

c i_ = P c M t 2 

rv! 



E — = ' 



Theorem 13.131 is proved. □ 



18 S. V. LOTOTSKY AND B. L. ROZOVSKII 

Remark 3.14. Taking M^w = c^u shows that, in general, bound A3. 33\) cannot be 
improved. When condition A3. 3(J\) does not hold, a bound similar to A3. 31}) can be 
established in a weighted space TZL 2 (¥;H), for example with r a = q a , where = 
l/(2 fc (l + Cfc)). For special operators a more delicate analysis might be possible; 
see, for example, [9]. 

If / and Mo are not deterministic, then the solution of (j3.4p might not satisfy 

E\\u(t)\\ 2 H <C(C A ,C M ,K,t) n E\\f(s)\\ 2 v ,ds + E\\u Q \f H 

even if all other conditions of Theorem 13.131 are fulfilled. An example can be con- 
structed similar to Example 9.7 in [18J: an interested reader can verify that the solu- 
tion of the equation u(t) = uq + J Q u(s)o£ ds, where £ is a standard Gaussian random 

variable and u = Va„ — satisfies Ew 2 (l) > — afa ■ For equations with 

n>0 n>l 

random input, one possibility is to use the spaces («S)_i i3 ; see (12. 6p . Examples of the 
corresponding results are Theorems 14.61 and 15. II below and Theorem 9.8 in [IB] . 

4. Stationary equations 

4.1. Definitions and Analysis. The objective of this section is to study stationary 
stochastic equation 

Au + 8(Mu) = f. (4.1) 

Definition 4.1. The solution of equation ( [^.1| ] with f G TZL 2 (¥;V') , is a random 
element u G 7ZL 2 (F;V) so that, for every ip satisfying ip G TZ~ 1 L 2 (E) and Dip G 
TZ~ l L 2 (¥; U), the equality 

((Au,<p)) + ((5(Mu),<p)) = ((f,<P)) ( 4 -2) 

holds in V . 

As with evolution equations, we fix an orthonormal basis it in U and use (I2.16P to 
rewrite ([4.1 p as 

Au + (Mm) oW = f, (4.3) 

where 

MuoW := ^M t tio( fc . (4.4) 

k>l 

Taking ip = £ a in (14.21) and using relation (12.171) we conclude, as in Theorem 13.51 that 
u = J2 a <=j u aia is a solution of equation (14. ip if and only if u a satisfies 

Au a + ^ \ f ®k M fc M Q _ £fe = f a (4.5) 

fc>l 

in the normal triple (V,H,V). This system of equation is lower-triangular and can 
be solved by induction on \a\. 

The following example illucidates the limitations on the "quality" of the solution of 
equation (14. ip . 
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Example 4.2. Consider equation 

u = l + uo£. (4.6) 

Similar to Example \3.b\ we write u = Yln>o u (n)Hn(Q/Vrd, where H n is Hermite 
polynomial of order n \2.4\) - Then ( |^.<5| ) implies «( n ) = J( n= o) + v^ u (n-i) or M( ) = 1, 
U(n) = Vn\, n > 1, or u = 1 + ^2 n>1 H n (0- Clearly, the series does not converge in 
L 2 (F) ; but does converge in (<5>)-i,,j for every q < (see Ii2.6\) ). As a result, even a 
simple stationary equation ^4-°] ) can be solved only in weighted spaces. 

Theorem 4.3. Consider equation ( [^.3] ) in which f G TZL 2 (¥; V) for some TZ. 

Assume that the deterministic equation AU = F is uniquely solvable in the normal 
triple (V, H, V), that is, for every F G V , there exists a unique solution U = A~ X F G 
V so that \\U\\v < C^H-Flly- Assume also that each is a bounded linear operator 
from V to V so that, for all v £ V 

HA^MfcuHv < C k \\v\\ v , (4.7) 

with Ck independent of v. 

Then there exists an operator TZ and a unique solution u G TZL 2 {¥\ V) of ^3. 4\ )- 

Proof. The argument is identical to the proof of Theorem 13.101 □ 

Remark 4.4. The assumption of the theorem about solvability of the deterministic 
equation holds if the operator A satisfies (Av, v) > n\\v\\y for every v G V, with k > 
independent of v. 

An analog of Theorem 13.111 exists if / is non-random. With no time variable, we 
introduce the following notation to write multiple integrals in the time-independent 
setting: 

where B is a bounded linear operator from V to V <S> U. 

Theorem 4.5. Under the assumptions of Theorem \4-3\ if f is non-random, then the 
following holds: 

(1) the coefficient u a , corresponding to the multi-index a with \a\ = n > 1 and 
the characteristic set K a = {k\, . . . , k n }, is given by 

7^7 $Z Bfe -w ' ' ' Bfc CT (D M (o)> ( 4 - 8 ) 

where 

• V n is the permutation group of the set (1, . . . , n); 

• B k = -A^M*; 

• U(o) = A" 1 /- 

(2) the operator TZ can be defined by the weights r a in the form 

a 00 

r Q = -±= where q a = J] q a k \ (4.9) 



'a! 
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where the numbers q k , k > 1 are chosen so that ^2 k>1 q k Cl < 1, and C k are 
defined in ffTfy - 
(3) With r a and q k defined by lj4-9\ ), 

X;? a «a€ B = ^ ) (A- 1 /), (4.10) 

\a\=n 

where B = — (giA -1 Mi, q 2 A~ 1 M. 2 , . . .), and 

Tin = A" 1 / + 4t ^(A-'fl (4.11) 



n>l 



Proof. While the proofs of Theorems 13.111 and 14.51 are similar, the complete absence 
of time makes equation (14. 3p different from either (13 .4p or anything considered in 
Accordingly, we present a complete proof. 

Define If / is deterministic, then v,(p) = A 1 f and, for \a\ > 1, 

Au a + E a k M k u a ^ ek = 0, 



fe>i 



or 



u. 



k>l k£K a 

where K a = {k\, . . . , k n } is the characteristic set of a and n = \a\. By induction on 



n. 



u. 



and (gSl) follows. 
Next, define 

= E n - °- 

\a\=n 

Let us first show that, for each n > 1, U n e L 2 (F; V). By (14. 8 p we have 



7/ II 2 < r : 

"ally ^ u - 



a! 



if 

fc>i 



By 



E^ a iW^^II*n! £ (^II(^) 

[a|=n |a|=n \ fc>l 



2a, 



(4.12) 



CiW v ,n\ 



E« 



< OO, 



,fe>l 



because of the selection of q k , and so C/„ G /^(F; V"). If the weights r a are defined by 
3, then 



2 II l|2 



E E a* ^ ^ii/iiv'E (e«) < °°> 

n>0 [ a |=n n>0 \fc>l / 
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because of the assumption ^2 k>1 C k q k < 1- 

Since (14.111) follows directly from (I4.10p . it remains to establish (I4.10p . that is, 

U n = <%(£/„_i), n > 1. (4.13) 

For n = 1 we have 

Ui = ^q k Ue k ik = ^B fc w (0) £ fc = frgiUo 



fe>i fe>i 



where the last equality follows from (I2.16p . More generally, for n > 1 we have by 
definition of C/ n that 

q a u a , if | ct | = n, 

0, otherwise. 



From the equation 



q a Au a + ^q k ^M k q a £k u a - £k = 
k>l 

we find 

{E V&k(lkB k q a ~ £k u a - £kl if |a| = n, 
k >i 
0, otherwise. 

= 2^ B fc (f/ n _i) Q _ £fe , 
fc>l 

and then (I4.13P follows from (I2.17p . Theorem 14.51 is proved. □ 

Here is another result about solvability of (14.31) . this time with random /. We use 
the space (<S) p , q , defined by the weights (12. 6p . 

Theorem 4.6. In addition to the assumptions of Theorem \4-3\ let Ca < 1 and 

C k < 1 for all k. If f G (S)-i^e(V') for some I > 1, then there exists a unique 
solution u G (<S)-i _^_4(V) of /jj4-3\ ) and 

IHIW- w (v)<C(^)||/||(5)_ 1 ,_ < tv)- ( 4 - 14 ) 

Proof. Denote by 11(^5 7), 7 G J, g 6 V, the solution of (14.31) with / Q = g/( a=7 ), 
and define « a = (a!)~ 1//2 M a . Clearly, u a (g,'j) = if |a| < I7I and so 

(4.15) 

It follows from (14.51) that 

u a+1 {Ul) =w Q (/ 7 ( 7 !)- 1/2 ;(0)). (4.16) 
Now we use (I4.12p to conclude that 

||W/r;7)llv<-fenll/lk'- (4-17) 
va'7! 
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Coming back to (14.151) with r 2 a = (a!)~ 1 (2N)^~^~ 4 ^ and using inequality (I2.2p we find: 

IH/-,;7)ll< 5) _„_«_ 4 <v> < C(t)(wt)- pp^ , 

where 




C(£) = 

( 12. 141) and ( 12.2ft imply C(£) < oo. Then (14. 14[) follows by the triangle inequality after 
summing over all 7 and using the Cauchy-Schwartz inequality. □ 

Remark 4.7. Example \4-$\ in which f G (<S)o,o and u G {S)-\ A , g < 0, shows that, 
while the results of Theorem \4-6] are not sharp, a bound of the type \\u\\(s) pq (y) < 
C||/||(5) p£ (V) is, in general, impossible if p > — 1 or q> I. 

4.2. Convergence to Stationary Solution. Let (V, H, V) be a normal triple of 
Hilbert spaces. Consider equation 

u{t) = (Au{t) + f{t)) + M k u{t) (4.18) 

where the operators A and do not depend on time, and assume that there exists 
an /* G 1ZL 2 (¥; H) such that lim^oo \\f(t) — f* \\-r,l 2 (W;H) = 0- The objective of this 
section is to study convergence, as t — > +00, of the solution of ( 14. 18f) to the solution 
u* of the stationary equation 

-Au* = f* + M k u*o£ k . (4.19) 

Theorem 4.8. Assume that 

(CI) Each Mfc is a bounded linear operator from H to H , and A is a bounded linear 
operator from V to V with the property 

(Av,v) + k\\v\\v <-c\\v\\ 2 H (4.20) 

for every v G V, with k > and c > both independent of v. 
(C2) / G TZL 2 i¥;H) and there exists an f* G 1ZL 2 {¥] H) such that 

lim^+oo \\f(t) - f*\\nL 2 (F;H)- 

Then, for every w G TZL 2 (¥; H), there exists an operator TZ so that 

(1) There exists a unique solution u G TZL 2 {¥\ V) of d4.18\ ), 

(2) There exists a unique solution u* G 1ZL 2 (W; V) of h4-19\ ), and 

(3) The following convergence holds: 

lim \\u(t) - u*\\ nL2{v , H) = 0. (4.21) 

Proof (1) Existence and uniqueness of the solution of (14. 18j) follow from Theorem 
EDI and Remark E2J 

(2) Existence and uniqueness of the solution of (14.191) follow from Theorem 14.31 and 
Remark IP1 

(3) The proof of (14.211) is based on the following result. 
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Lemma 4.9. Assume that the operator A satisfies U^Wty and F = F(t) is a de 



terministic function such that lim 4 ^ +00 = 0. Then, for every U G H, the 

solution U — U(t) of the equation U(t) — Uq + J * AU{s)ds + J" * F(s)ds satisfies 
lim^ +00 \\U(t)\\ H = 0. 

Proof. If $ = $ 4 is the semi-group generated by the operator A (which exists because 
of flCTD ), then 

U(t) = ® t U + / ® t -sF(s)ds. 
Jo 

Condition (14.201) implies ||$tC/ ||i? < e' ct \\U \\ H , and then 

\\U(t)\\ H < e- ct \\U \\ H + t e- c{t -^\\F(s)\\ H ds. 



o 



The convergence of \\U(t) \\h to zero now follows from the Toeplitz lemma (see Lemma 
A. 2 in Appendix). Lemma [4.91 is proved. □ 

To complete the proof of Theorem I4.8[ we define v a (t) = u a (t) — u* a and note that 

k 

By Theorem 14. 3^ u* a G V and so v a (0) G H for every a G J . By Lemma 14.9} 
lim^ +00 ||f(o)(t)||j? = 0. Using induction on \a\ and the inequality ||M fc t> Q _ £fe (t)||# < 
Ck\\v a - ek (t)\\H, we conclude that lim^+oo ||f a (i)llH — for every a G J. Since 
v a G C((0,T); H) for every T, it follows that sup t>0 ||i> a (i)||# < oo. Define the 
operator 71 on L 2 (F) so that lZ£, a = r a £ a , where 

(2N)~ Q 



l + sup\\v a {t)\\ H 

t>o 

Then (14.211) follows by the dominated convergence theorem. 

Theorem 14.81 is proved. □ 

5. Bilinear parabolic and elliptic SPDEs 

Let G be a smooth bounded domain in IR d and {hk, k > 1}, an orthonormal basis in 
L 2 (G). We assume that 

sup \h k (x)\ < Cfc, k > 1. (5.1) 
A space white noise on L 2 (G) is a formal series 

W(x) = Y,hk(x)£ k , (5.2) 

k>l 

where k > 1, are independent standard Gaussian random variables. 
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5.1. Dirichlet Problem for parabolic SPDE of the Second Order. Consider 
the following equation: 

u t (t,x) = aij(x)DiDjii(t,x) + bi(x)Diu(t,x) + c(x)u(t,x) + f(t,x) 

(5.3) 

+ (ai(x)Diu(t, x) + u(x)u(t, x) + g(t, x)) o W(x), < t < T, x G G, 

with zero boundary conditions and some initial condition u(0,x) = Uq(x); the func- 
tions Oy, bi, c, /, <Tj, v, g, and Uq are non-random. In (I5.3P and in similar expressions 
below we assume summation over the repeated indices. Let (V, H, V) be the normal 

triple with V = H\(G), H = L 2 (G), V = H 2 1 (G). In view of ([ED, equation Q is 
a particular case of equation (13 .4p so that 

Am = a i j{x)D i DjU + bi(x)DiU + c(x)u, M fc u = (ai(x)DiU + u(x)u)hk(x), (5.4) 

and f(t, x) + g(t, x) o W{x) is the free term. 

We make the following assumptions about the coefficients: 

Dl The functions Lipschitz continuous in the closure G of G, and the 

functions bi, c, Oi, v are bounded and measurable in G. 
D2 There exist positive numbers A\,A 2 so that 1 2/ 1 2 < a-ij{x)yiyj < ^bM 2 for 

all x G G and y G M d . 

Given a T > 0, recall the notation V = L 2 {{Q,T); V) and similarly for 7i and V (see 

(SU)). 

Theorem 5.1. Under the assumptions Dl and D2, if f G V , g G 7Y ; m G 

t/ien there exists an i > 1 and a number C > 0, both independent of Uq, f,g, so that 
u G 1ZL 2 (¥; V) and 

IMkL 2 (F;V) < C • (||«o||h + \\f\\v + \\g\\n), (5-5) 
where the operator TZ is defined by the weights 

r 2 a = c - 2a (\a\\)- 1 (2Ny 2ea (5.6) 
and c a = rifc c fc fe ; with Ck from ( 15. the number i in general depends on T . 

Proof. We derive the result from Theorem 13. Ill Consider the deterministic equation 
U(t) = AU(t) + F. Assumptions Dl and D2 imply that there exists a unique solution 
of this equation in the normal triple (V, H, V), and the solution satisfies 

sup \\U(t)\\ H + \\U\\ V < C ■ (\\U(0)\\h + \\F\\v>), (5.7) 

0<t<T 

where the number C depends on T and the operator A. Moreover, 05. ip implies 
that (I3.15P holds with Ck = C Ck for some positive number C independent of k, but 
possibly depending on T. 

To proceed, let us assume first that g = 0. Then the statement of the theorem follows 
directly from Theorem 13 .111 if we take in (I3.18P qt = c^" 1 (2A;)~ £ with sufficiently large 

L 
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It now remains to consider the case g 7^ and / = Uq = 0. Even though g is non- 
random, g^k is, and therefore a direct application of Theorem 13.111 is not possible. 
Instead, let us look more closely at the corresponding equations for u a . For a = (0), 

U( )(t) = / Aii(o)(s)ds, 
Jo 

which implies u^(t) = for all t. For a = 

u £k (t) = / Au £k (s)ds + h k \ g(s)ds, 
Jo Jo 

or 



u £k (t) = / ^ t -sh k g(s)ds, 
Jo 



so that 

IKJv < CoCjfcll^Hw. (5.8) 

If \a\ > 1, then 



U a (t) = / Attq(s)d.S + ^ y^Mfcttg,-^, 

- 70 fc>i 



which is the same as ( 13.22}) . In particular, if \a\ = 2 and {i, j} is the characteristic 
set of a, then 

i rt 



(*) = -=/ (M<u ei (s) + M> e4 (s)) ds. 



More generally, by analogy with (I3.25p . if |ct| = n > 2 and {fci, . . . , k n } is the charac- 
teristic set of a, then 



1 . p f Sn f S3 

Ua(t) = —j= J2 / • • • / $ *-^ M fe CT( „) • ■ • ®s 3 -s 2 M ka{2) U £a(1) (s2)dS2 . . . rfs n . 

By the triangle inequality and (15. 8p . 



«a V < 7= IMIw, 



Va! 

and then (15.51) follows from (12.21) if I is sufficiently large. 

This completes the proof of Theorem 15.11 □ 
Theorem 5.2. In addition to Dl and D2, assume that 



(1) o~i = for all i; 



(2) the operator A in G with zero boundary conditions satisfies ( 4-20 ). 

If there exist functions f* and g* from H so that 

lim (\\f(t)-r\\H+\\g(t)-g*\\ H ) = J (5.9) 

t— >+oo 

then the solution u of equation ( 15.31) satisfies 

lim \\u(t) - u*\\tzl 2 (V;H) = 0, (5.10) 

£ — >+co 
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where the operator 1Z is defined by the weights ( f5. bty and u* is the solution of the 
stationary equation 

aij(x)DiDjU*(x) + bi(x)DiU*(x) + c(x)u*(x) + f*(x) 

+ (u(x)u*(x)+g*(x))oW(x) = 0, x G G; u\ dG = 0. < ' ) '" ) 

Proof. This follows from Theorem 14.81 □ 

Remark 5.3. The operator A satisfies l[4-20\ ) if, for example, each a%j is twice con- 
tinuously differentiable in G, each bi continuously differentiable in G, and 

inf c(x) - sup(|£VDj-ay(a:)| + |A&i(^)|) > £ > 0; (5.12) 

^eG x( zg 

this is verified directly using integration by parts. 

5.2. Elliptic SPDEs of the full second order. Consider the following Dirichlet 
problem: 

-A (x) Dju {x)j + 
Di (aij (x) Dj (u (x)) ) o W (x) =f(x), xeG, ( 5 - 13 ) 
u\dG = 0, 

where W is the space white noise (15. 2p . Assume that the functions a^, a^, f, and g 
are non-random. Recall that according to our summation convention, in (15 . 13[) and 
in similar expressions below we assume summation over the repeated indices. 

We make the following assumptions: 

El: The functions ) and measurable and bounded in 

the closure G of G. 

E2: There exist positive numbers A\,A% so that v4i|y| 2 < a-ij{x)yiyj < A 2 \y\ 2 

for all x G G and y G W 1 . 
E3: The functions h k in (15.21) are bounded and Lipschitz continuous. 

Clearly, equation (15. 13j) is a particular case of equation (14. 3 \ with 

Au(x) := -Di (dij (x) DjU (x) J (5.14) 

and 

Wl k u{x) := h k (x) Di^cTij (x) DjU (x) j. (5.15) 
Assumptions El and E3 imply that each is a bounded linear operator from 

o 

H 2 1 (G) to H^iG). Moreover, it is a standard fact that under the assumptions El 
and E2 the operator A is an isomorphism from V onto V (see e.g. [17]). Therefore, 
for every k there exists a positive number C k such that 

\\A- l M k v\\ v <C k \\v\\ v , veV. (5.16) 
Theorem 5.4. Under the assumptions El and E2, if f G H^iG), then there exists 

o 

a unique solution of the Dirichlet problem $5.1 Sty u G 7^X2 (F; H\(G)) such that 

\\u\\ . <C- 11/11 H -i (rv (5.17) 
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The weights r a can be taken in the form 



q 



a 



OO 



where q a = JJg^ fc , (5.18) 



and the numbers q^, k > 1 are chosen so that J\>i ^klk < 1? with Ck from 115. 16]) . 

Proof. This follows from Theorem 14. 5[ □ 
Remark 5.5. With an appropriate change of the boundary conditions, and with extra 



regularity of the basis functions hk, the results of Theorem 5.4 can be extended to 
stochastic elliptic equations of order 2m. The corresponding operators are 

Au = (-1)™ A x • • • D im ("„..,..,,..„ Or) D n ■ ■ ■ D jm u (x) ) (5.19) 

and 

M k u = h k (x) D h ■ - ■ D im (o h „ Amh „, jm (x) D h ■ ■ ■ D jm u (x) ) . (5.20) 

Since G is a smooth bounded domain, regularity of h k is not a problem: we can take 
h k as the eigenf unctions of the Dirichlet Laplacian in G. 
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Appendix. 

A FACTORIAL INEQUALITY. 

Lemma A.l. For every multi-index a G J , 

\a\\ < a\(2N) 2a . (Al) 

Proof. Recall that, for a — (a 1: . . . , a k ) G J , 

k k k 

\a\ = J>, a\ = l[a f \, N a = \\t a * 

i=\ e=i i=\ 

It is therefore clear that, if \a\ = n, then it is enough to establish (Al) for a with 
«fc = for k > n + 1, because a shift of a multi-index entry to the right increases 
the right-hand side of (Al) but does not change the left-had side. For example, if 
a = (1, 3, 2, 0, . . .) and f3 = (1, 3, 0, 2, 0, . . .), then \a\ = a\ = but N Q < N beta . 
Then 

- I 2 2 v? ^ a] N 2a ' 

x ' a 1 +...+a n =n 



STOCHASTIC EQUATIONS WITH SPATIAL NOISE 29 

where the equality follows by the multinomial formula. Since all the term in the sum 
are non- negative, we get (A2). 

□ 

The proof shows that inequality (A2) can be improved by observing that J2k>i k~ 2 — 
7r 2 /6 < 2. One can also consider J2 k>1 k~ q for some 1 < q < 2. 

A VERSION OF THE TOEPLITZ LEMMA. 

Lemma A. 2. Assume that f = f(t) is an integrable function and\im t _ >+OQ \f(t)\ = 0. 
Then, for every c > 0, lim t ^. +00 J* * e~ c ^~ s ^f(s)ds = 0. 

Proof. Given e > 0, choose T so that < e for all t > T. Then 

\Joe-< t -^f(s)ds\ < e~ ct J T e cs \f(s)\ds + e fee'^-^ds. Passing to the limit as 
t — > +oo, we find lim t _ >+00 | J* e~ c ^~ s ^f(s)ds\ < e/c, which completes the proof. 

□ 
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